Non-gaussianity in axion N-flation models: detailed predictions and mass spectra 
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We have recently shown [lj that multi-field axion N-flation can lead to observable non-gaussianity 
in much of its parameter range, with the assisted inflation mechanism ensuring that the density 
perturbations are sufficiently close to scale invariance. In this paper we extend our analysis in several 
directions. In the case of equal-mass axions, we compute the probability distributions of observables 
and their correlations across the parameter space. We examine the case of unequal masses, and show 
that the mass spectrum must be very densely packed if the model is to remain in agreement with 
observations. The model makes specific testable predictions for all major perturbative observables, 
namely the spectral index, tensor-to-scalar ratio, bispectrum, and trispectrum. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

In a recent paper [l[ we identified a new mechanism 
for generating an observably large non-gaussianities dur- 
ing inflation. The effect is due to diverging trajectories 
near a maximum of the potential, where — in a single- 
field model — the field would have a large effective mass 
in Hubble units. Such large masses are excluded in the 
single-field case because the scalar spectral index is too 
far from unity, but in the multi-field case the assisted 
inflation phenomenon Q generates a spectrum that can 
be compatible with current observational constraints Q 
without simultaneous suppression of the dimensionless 
bi- or tri-spectra. As a specific example, we imple- 
mented this mechanism in the multi-field axion N-flation 
model Q. An alternative implementation of the mecha- 
nism in a hybrid inflation context was recently given by 
Mulryne et al. Q. 

The purpose of the present paper is to provide a 
more detailed phenomenological description of the axion 
N-flation model. In the case where all fields have equal 
mass, which was assumed throughout Ref. [l|, we provide 
a comprehensive analysis of all key observables, analyz- 
ing the probability distributions inherited from random- 
ness in the initial conditions, their dependence on model 
parameters, and their intercorrelations. We extend the 
analysis to the case where the fields are distributed with 
a spectrum of masses, and obtain tight constraints on 
the packing of the mass fraction necessary to maintain 
agreement with observations. 



II. THE AXION N-FLATION MODEL 

The axion N-flation model is based on a set of N{ un- 
coupled fields, labelled <fii, each with a potential Q 
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the fields but we will not consider these. The mass 
of each field in the minimum of the potential satisfies 
rrii = 27rA|//j, and the angular field variables a.% lie in 
the range (— 7r, +7r] . Without loss of generality we will set 
initial conditions with all on positive. If only a single field 
is present this model is known as natural inflation . 

One motivation for N-flation was to avoid the require- 
ment for super-Planckian field values Q, which are in- 
voked in many single-field models. If one literally im- 
poses \<f>\ < Mp (where Mp = (87rG) -1 / 2 is the reduced 
Planck mass) this requires fi < 2Mp for each i. How- 
ever, it would be reasonable to regard this condition as 
a guideline rather than mandatory. 



A. The amount of inflation 

Any inflationary model must provide sufficient 
e-foldings to resolve the classical cosmological problems. 
For a given set of initial angles a* one finds 
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where on = 2ir(j)i/fi and fi is the i th axion decay con- 
stant. More generally, couplings may exist between 



where in the second line we have ignored a small cor- 
rection from the location of the end of inflation. The 
sum is dominated by fields whose initial angle is close 
to 7r. However, for any reasonable distribution of a*, 
the logarithm means that many fields must cooperate to 
yield sufficient e-foldings — unless / is extremely large in 
Planck units, as in natural inflation. As remarked above, 
this conflicts with the goal of maintaining sub-Planckian 
field excursions. Our successful models typically feature 
hundreds or thousands of fields. 

If the initial conditions are taken to be distributed uni- 
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formly in angle, one can average to find 
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which confirms the requirement for many fields. 

Having obtained sufficient inflation, we must identify 
the epoch at which observable perturbations were gener- 
ated. This requires knowledge of the entire history of the 
Universe, including reheating, and is therefore subject to 
some uncertainty. We follow Ref. to obtain the time 
at which the present horizon scale k = aoH crossed out- 
side the horizon, expressed in e-foldings before the end 
of inflation: 
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The last term might typically be —5 0j- Here 'hor', 'reh' 
and 'end' denote values at horizon crossing and at the end 
of reheating and inflation, respectively. To evaluate these 
terms, we must ensure that models are normalized to 
reproduce the correct amplitude of density perturbations. 

For a single field, the appropriate JVhor depends how 
close the initial position lies to the maximum; we find 



iVhor — 54 (near the maximum) , 

~ 59 (away from the maximum) . 



(5) 



Near the hilltop, we require fewer e-folds because the 
flatness of the potential implies we require smaller H*. 
Therefore, the duration of the radiation era which follows 
inflation is shorter. 

We aim to verify that iVhor does not shift significantly 
when Nf 3> 1. Taking Nf of order 10 3 , we find 
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We have assumed / = Mp but the result is relatively 
insensitive to this choice. We conclude there is no signif- 
icant change to Nh or . 

Finally, we note that the constraints we impose are not 
evaluated at the present-day horizon scale 
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but rather at fc* = 0.002 Mpc - . This scale is inside 
the present horizon, with iV* ~ iVhor — 2. Overall, we 
conclude that in the multi-field case it remains a reason- 
able hypothesis that the pivot scale crossed the horizon 
around 50 to 60 e-foldings before the end of inflation. In 
what follows we will consider only these two values. 



B. Perturbations 

We calculate observables using the SN formula [§], 
which measures fluctuations in the total e-foldings of ex- 
pansion owing to field perturbations. We define e-like 



slow-roll parameters for each field, 



(8) 



where a prime denotes a derivative with respect to 0j, 
and no summation over i is implied. The global slow- 
roll parameter e = —H/H 2 can be written as a weighted 
sum e ~ ^ZiiVi/V) 2 ^, in which each field contributes 
according to its share of the total energy density. We 
must have e < 1 during inflation. 

We work in the horizon-crossing approximation, in 
which the dominant contribution to each observable is 
assumed to arise from fluctuations present only a few 
e-folds after horizon exit of the wavenumber under dis- 
cussion. After smoothing the universe on a superhori- 
zon scale somewhat smaller than any scale of interest, 
the horizon-crossing approximation becomes valid when- 
ever the ensemble of trajectories followed by smoothed 
patches of the universe approaches an attractor. The va- 
lidity of the horizon-crossing approximation was recently 
discussed by Elliston et al. j9|, the expectation being 
that it is a good approximation in our case. Numeri- 
cal calculations supporting this conclusion were reported 
in Ref. 0. 

The observables of interest, defined in the conventional 
way Hd are the scalar spectral index n, the tensor-to- 
scalar ratio r, the bispectrum /nl, and the trispectrum 
parameters tnl and c/nl- They are given by 
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where TV.;, N ti j and N t ij k are respectively the first, sec- 
ond and third derivatives of TV with respect to field val- 
ues at time *, corresponding to evaluation at the pivot 
scale determined as in Eq. $Z§. In writing Eqs. (JT^-fHU), 
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any intrinsic non-gaussianity among the field perturba- 
tions at horizon crossing has been neglected. This is a 
good approximation whenever the bi- and tri-spectrum 
parameters are large enough to be observable p^4l5| . 
Our sign convention for /nl matches that used in WMAP 
papers 0] , and the non-gaussianity is predicted to be of 
local type. The observed amplitude of perturbations is 
obtained by adjusting the A, to give an appropriate value 
of IT*. 

Under a quadratic approximation to each potential, it 
can be shown that Eqs. (fTTj) and (1X21) recover their single- 
field values of order ~ 1/AT* making /nl unde- 
tectably small. The spectral index can be shown to be 
less than its single- field value 1 — 2/iV* [I?} with equality 
only in the equal-mass case. Its value for a given choice 
of parameters must be computed numerically 18 1 . How- 



ever these results change whenever the initial conditions 
populate the hilltop region. 



III. THE EQUAL-MASS CASE 

In Ref. [H we considered only the case where all fields 
have the same potential. In this section we make the 
same assumption, but carry out a much more detailed 
analysis of the phenomenology. The scale A = A, is fixed 
from the observed amplitude of Vq, leaving f = fa and 
Nf as adjustable parameters. 

The initial conditions are drawn randomly from a uni- 
form distribution of angles en , with several realizations to 
explore the probabilistic spread. This choice seems plau- 
sible in light of the approximate circular symmetry of ax- 
ion potentials. One could of course envisage other prob- 
ability distributions, but to obtain a successful model 
there must be some reasonable probability of populat- 
ing the hilltop region. Since those fields dominate the 
statistics of the density perturbations, any probability 
distribution which is approximately flat near the hilltop 
can be expected to yield similar results. 



A. Analytic approximations 

The ti approach zero for fields close to the hilltop, so 
each summation in Eqs. ([9])-(|14[l is dominated by those 
fields with the smallest £j. Suppose some number N of 
such fields have roughly comparable , of order e. 

Quantum diffusion. Near the hilltop, the parameters 
Ci ~ e are small and the classical motion of each axion 
becomes small. In this region, two new effects emerge. 
First, for sufficiently small e the classical motion of indi- 
vidual fields can be dominated by quantum fluctuations, 
but typically this is not of concern unless the fields in- 
volved contribute non-negligibly to the energy density. 
In this case one can expect density fluctuations of order 
unity, leading to a phase of 'topological inflation' [T^j. 
Second, the hilltop is a singularity of the e-folding his- 
tory, N, as a function of the initial field values. If these 



initial values are chosen too close to the singularity then 
the Taylor expansion used to obtain Eqs. ([9])— (|14JI be- 
comes unreliable. 

We consider the constraints in turn, beginning with the 
issue of singularities in N. These emerge from Eq. @ 
in the limit a* — > tt. By repeated differentiation, we 
conclude that the Taylor expansion is trustworthy unless 
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where the final approximate equality applies in a model 
for which / ~ Mp, and the field fluctuation 5(f) should be 
estimated at the time of horizon exit, which we continue 
to label *. Assuming N axions dominate the spectrum 
with comparable 5* ~ <5* and /, ~ /, the observed am- 
plitude of density fluctuations Vq — 2 x 1CU 9 requires 
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For N = O(10), Eqs. (fT5 ]) -(fTp ]l imply that a breakdown 
of the perturbative SN formula cannot occur unless at 
least one fa is a few orders of magnitude less than the 
Planck scale. For JV > 1 a more extreme tuning of some 
fi is required. As explained above, after drawing initial 
conditions a* within our numerical simulations, we ad- 
just the Hubble scale to satisfy Eq. (| 1C>[) by a suitable 
normalization of the scales A,;. 

Under normal circumstances the power spectrum 
is monotonically increasing with time and therefore 
Eq. (fl6|) guarantees that the adiabatic trajectory is stable 
to quantum fluctuations. But it should also be checked 
that when the final axion rolls to its minimum, forcing the 
correlation functions of £ to their horizon-crossing values 
Eqs. ©-(HU, its fluctuations are not large enough to ini- 
tiate an unwanted phase of topological inflation. There- 
fore we require 
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where in the final step we have taken / ~ Mp and, be- 
cause the final axion field dominates the potential by def- 
inition, we have estimated H m \\ ~ K 2 /Mp. 

Our numerical simulations do not take quantum diffu- 
sion into account, so we will usually wish to impose the 
stronger requirement that diffusion does not occur for any 
field. This also avoids the possibility that the final axion 
field diffuses to sufficiently small values that Eq. (|T7|) is 
violated. It is sufficient to demand that Eq. ([To]) bounds S 
away from the quantum diffusion regime. Taking / ~ Mp 
this requires 



(18) 



Once the spectrum has been correctly normalized, 
Eq. <118p can be interpreted as a bound on the number of 
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axions, TV c ff which contribute an energy density of order 
A 4 , 



N cS < N 1 '^- 1 ' 2 



(19) 



Note that -/V e ff > TV, since an axion which contributes 
an energy density of order A 4 will not contribute to 7\ 
unless its contribution is enhanced by proximity to the 
hilltop. A similar discussion was given by Huang [2(ij] . 

Observable quantities. We now proceed to study the 
various observable quantities, each of which has a dif- 
ferent scaling with TV. The spectrum, 7\, scales like TV 
copies of a single-field model with slow-roll parameter e, 
whereas r is reduced by a factor TV compared to its value 
in the same single-field model. The spectral index can be 
written 



n — 1 ~ — 2e* — 
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and is independent of TV. Instead, the summation in the 
denominator receives contributions from all fields. As- 
suming / is not too different from Mp, the spectral in- 
dex becomes close to — 2e» when this sum is of order 10 3 . 
This is the familiar assisted-inflation mechanism. In con- 
trast, the bispectrum amplitude /nl has the approximate 
behaviour 
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and is independent of e if the dominant fields are suf- 
ficiently close to the hilltop. It is the different scalings 
of n — 1 and /nl with TV which makes the scenario vi- 
able: the N-flation mechanism lifts the sing le-field con- 
sistency condition /nl ~ — (5/12)(n — 1) Qjl], which pre- 
vents single-field models generating large non-gaussianity 
without violating observational bounds on n. 

A similar analysis applies to the trispectrum, for which 
it is conventional to parameterize the amplitude of a 
local-type trispectrum using the parameters tnl and <7nl 
of Eqs. pll-rtlgll [U, 
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Where the summations in Eqs. (|9|)- (fT2[) are dominated 
by a single field, this formula shows that the non-gaussian 
parameters can become rather large, scaling as powers of 
(M P //) 2 . For / = M P , we find / NL < 16.4. A non- 
gaussian fraction of this magnitude should be visible to 
the Planck satellite. The same parameter choice yields 
tnl ;$ 390 and #nl 360. Such a small #nl is unlikely 
to be observable, although there is some hope that tnl 
of this order could be detected with a future microwave 
background polarization satellite [22j . 
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FIG. 1: Predictions in the n-r plane, averaged over realiza- 
tions, for various values of / between 0.4Mp and 2Mp and of 
TVf between 464 and 10,000, all giving sufficient inflation. The 
black (left) cluster of points takes TV, = 50 and the red (right) 
cluster TV, = 60. The quadratic expansion predicts r = 8 /TV, , 
far off the top of this plot. The region right of the line is 
within the WMAP7+BAO+#o 95% confidence contour [1. 



For smaller / the density perturbation becomes in- 
creasingly non-gaussian. It is even possible to achieve 
/nl ~ 100 for / ~ 0.4M P , although then TVf must be 
very large to gain sufficient e-foldings and some tension 
with Eq. (fT9|) may emerge. At this /, the trispectrum 
parameters may become as large as tnl ~ 1-5 x 10 and 
5nl ~ 1.4 x 10 4 . Eq. (f2"2")) shows that when a few hill- 
top fields dominate the density perturbation, the axion 
N-flation model reproduces the single-field relation be- 
tween /n t, an d tnl first pointed out by Suyama and Ya- 
maguchi 



B. Numerical calculations 

To fully explore the model space requires numerical 
calculations of the evolution, which we carry out using 
an extension of the code developed in Ref. [18| . For each 
choice of model parameters, a set of runs is required to 
explore the uncertainty induced by the random initial 
conditions. 

In Fig. [T] we show model predictions in the n-r plane, 
averaged over several realizations of the initial condi- 
tions. We see n and r are only weakly dependent on the 
model parameters (though there is significant dispersion 
amongst realizations, not shown here), with the choice of 
TV, being the principal determinant of n. The models are 
compatible with current observational constraints in the 
n-r plane. 

Turning to the non-gaussianity, Fig. [5] shows /nl as 
a function of TVf for / = Mp, with ten realizations at 
each TVf. This clearly shows the expected maximum, 
which is nearly saturated in cases where a single field 
dominates the summations. In cases where several fields 
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FIG. 2: Predicted non-gaussianity, |/nl, for / = Mp and FIG. 4: The distribution of the tensor-to-scalar ratio, r, due 
TV, = 50. The error bars are on the mean over realizations to initial condition randomness, for / = Mp and N* = 50. 
(not the standard deviation). Here the maximum achievable 
value of I /nl is 2tt 2 ~ 20, almost saturated in some real- 
izations. The significant spread is due to initial condition 
randomness with typical mean values being around half the 
maximum achievable value, and no discernible trend with Nf. 
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FIG. 3: The distribution of non-gaussianity, (6/5)/nl due to 
initial condition randomness, for / = Mp and TV, = 50. 



contribute significantly to the sums in Eqs. (|9|)- (|T2l) . the 
non-gaussian fraction is reduced. 

It is clear from Fig. [2] that the non-gaussianity has a 
large variance between different realizations of the ini- 
tial conditions. To study this in more detail, Fig. [3] 
shows the distribution of /nl for / = Mp and AT* = 50 
with N{ = 2150, now for 100 initial condition realiza- 
tions. Values near the maximum, corresponding to the 
non-gaussianity signal being dominated by a single field, 
occur about 25% of the time, and then there is a broad 
peak at smaller values indicating an effective number of 
contributing fields around two or more. The broad dis- 
tribution implies that a measurement of /nl alone could 
not accurately constrain model parameters. 
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FIG. 5: A scatter plot between the non-gaussianity and the 
tensor-to-scalar ratio, /nl^t due to initial condition random- 
ness, for / = Mp and N, = 50. 



Fig. [T] shows that we always obtain r values much 
smaller than the 8/./V* predicted in single- field models. 
Again there is significant dispersion from initial condi- 
tions, shown in Fig. [U with the actual values making 
any observation challenging in the extreme. Note that 
Fig. H] is with the same conditions as Fig. [3] 

Moreover, /nl and r are strongly correlated, as shown 
in Fig. [5l with the larger values of /nl corresponding to 
smaller ones of r. The interpretation is that large /nl 
requires one of the fields to be very near the maximum, 
where the flat potential forces down the normalization A^ 
which then takes r down as well. 

In Fig. we show the predicted non-gaussianity as a 
function of /, for a range of choices of Nf. Each point 
shown is the average of five or more realizations for an 
f-N{ pair. We see a strong trend with /, well represented 
by Eq. (JUJ) with TV ~ 2. The different N { arc scattered 
by randomness in the initial conditions rather than an 
identifiable trend. 
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FIG. 6: The predicted non-gaussianity as a function of /, for 
various N{. 
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FIG. 8: A scatterplot of the bispectrum /nl versus the 
trispectra due to initial condition randomness, for / = Mp 
and N* = 50. Black crosses (x) and blue circles (o) are for tnl 
and (54/25)<7nl respectively. The red line shows the single- 
field relation between /nl and tnl from Ref. [2^ ]. 
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(see also Ref. [2J] for a more general derivation) . We see 
that our models do indeed satisfy this inequality on a 
case-by-case basis, and approach equality in the limit of 
the highest achievable non-gaussianity. 
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IV. UNEQUAL MASSES 
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FIG. 7: The predicted trispectrum tnl as a function of /, for 
various Nf. 



Fig. [7J shows the prediction of the trispectrum of non- 
gaussianity tnl, as a function of /, for a range of choices 
of JVf as well as ones in /nl- Again each point shown is 
the average of five or more realizations for an f—Nf pair, 
and again we see the strong trend with /, well matched 
by Eq. $Z2§ with N ~ 2. We have not plotted the cor- 
responding figure for (54/25) <7nl because it is so similar 
to tnl that it would be almost identical to Fig. [7J This 
also follows the trend of Eq. (|2"5j) . 

To study the relation between bispectrum and trispec- 
trum, Fig. [5] shows the predicted trispectrum as a func- 
tion of / NL for / = M P and = 50 with N { = 2150. 
The relation shows much less scatter than do the individ- 
ual quantities, due to their origin in common dynamics. 
For large /nl, the scatter is at its smallest, because one 
field needs to dominate in this case and this field gen- 
erates each non-gaussiani ty p arameter in the same way. 
Suyama and Yamaguchi [23[ demonstrated that single- 
field models satisfy tnl = [(6/5)/nl] 2 , and that in more 
general models this expression gives a lower bound to tnl 



We consider now the unequal-mass cases. We take the 
mass spectrum as exponentially distributed: 



m exp 
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for i = 1,2, 



N { , (24) 



where m is the smallest mass. We studied this mass 
spectrum in the quadratic potential case in Ref. fl8| . Our 
main objective in this section is to constrain a, which 
governs how tightly the mass spectrum is packed. 

Choosing the mass spectrum does not fix the model, 
because the potentials depend on two parameters which 
combine to give the mass. We consider two extreme pos- 
sibilities. One is varying the amplitude of the potential 
of each field, A,, by giving the same value of the con- 
stant decay fi = / to each field, and the other is varying 
the decay constant fi while keeping the same amplitude 
Aj = A. In the former case the potentials all have the 
same period but different amplitudes, and in the latter 
the same amplitudes and different periods. 



A. Varying Aj 



Varying A, with fixed fi = f requires 
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FIG. 9: Predictions in the n-r plane for a = 500, 2000, 10000 
from top to bottom with various / up to a maximum / = Mp. 
The black (left) cluster of crosses takes N» = 50 and the red 
(right) cluster of circles N* = 60. Observationally allowed 
models lie to the right of the line. 



If we adopt this in Eqs. (|9]) to (fT4|) . then we see that the 
effect from the different amplitude acts only on the spec- 
tral index and not on r or the non-gaussianity parameters 
(/nl,Tnl, and g NL ). 

In Fig. 121 we show the predictions for n and r for a = 
500, 2000, and 10000 with / = M P . For a = 500, the 



FIG. 10: Locations in the a-f plane for N* = 50 (top) and 
N* = 60 (bottom) where observable predictions for the aver- 
age value n axe within the limits (circles, o) or outside them 
(crosses, x). 



numerical runs have been done with Nf in the range 1470 
to 6810, for a = 2000 with N f from 1470 to 10000, and 
for a = 10000 case with N f from 2150 to 10000. Roughly 
speaking, for a given a the obscrvablcs are independent 
of Nf. The spectral index n depends significantly on a, 
but we found the dependence on / is weaker. 

Small a forces the spectral index outside of its allowed 
region in the majority of cases, leading to a lower limit 
on a. Looking at various / and a, we can map out the 
allowed parameters. Figure [TU] samples the a-f plane 
to determine where the mean value of the spectral in- 
dex n exceeds the 95% observational limit n = 0.93, for 
cases with Nf in the range 464 to 10000. Even though 
each point is an average over ten initial condition real- 
izations (as many as we could reasonably run), there is 
still residual noise meaning there is not a perfect parti- 
tion of the parameter space into allowed and disallowed 
regions. Nevertheless, the trend is clear; small a is dis- 
favoured, while for large enough a the equal-mass limit is 
effectively attained which we already know to be viable. 
Match with data is achieved more comfortably for the 
larger iV» choice. 
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FIG. 11: Predictions in the n-r plane for a = 2000 (top) and 
10000 (bottom). The black crosses (x) take JV» = 50 and the 
red circles (o) N* = 60. 



The main result from varying Aj is that the mass spec- 
trum has to be tightly packed, and then the results from 
the equal-mass case are recovered. In all these cases there 
is no difference in the non-gaussianitics; with the same / 
the /nl does not change. 



B. Varying f t 



The second case fixes A» = A, implying 



fi 



2ttA 2 



(26) 



FIG. 12: Predictions in the a-A plane for N* = 50 (top) 
and N, = 60 (bottom), for several choices of Nf. Each point 
shown is the average of ten realizations for an a—Nf pair. Cir- 
cles (o) denote that the observable predictions for the spectral 
index are within the limits, and crosses (x) that the values are 
excluded. 



/ whose normalization depends on the particular ini- 
tial condition realization. Nevertheless, the same general 
trend is apparent that small a is disfavoured, with no 
working models found for a = 2000 regardless of AT* . For 
high enough a the models are always allowed, and in the 
intermediate regime their validity is a matter for detailed 
individual analysis. Once we restrict to observationally- 
allowed models, the non-gaussianity is essentially that 
of the equal-mass case though for a smaller 'effective' / 
value given the spectrum of fi values. 



where we constrain that the largest fi is always Mp . Nu- 
merical calculations were done with Nf in the range 2150 
to 10000 and a from 2000 to 10000. There are not enough 
e-foldings in cases with a < 2000. The results in the n-r 
plane, and the regions of parameter space where viable 
values of those are achieved, are shown in Figs. QT] and 

The parameter space limits are less clear than in the 
previous case, as A is a less fundamental parameter than 



V. CONCLUSIONS 

We have carried out a detailed study of the phe- 
nomenology of the axion N-flation model, extending our 
previous work in several directions. This includes exten- 
sion of non-gaussianity calculations to the trispectrum, 
and an analysis of the correlations between different ob- 
servables induced by the initial condition realizations. 
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When a spectrum of unequal masses is considered, we 
find that the spectrum must be extremely tightly packed 
if the spectral index is to stay in agreement with obser- 
vations. This echoes the result we found for quadratic 
potentials in Ref. [l8| . Once this condition is obeyed, we 
find that the predictions for other observable quantities 
essentially match those of the equal-mass case, i.e. the 
mass spectrum does not introduce any new phenomenol- 
ogy 

Nevertheless, the model is highly predictive in terms 
of all the major perturbation observables, with the spec- 
tral index already close to the observational lower limit 
and the non-gaussianity detectable across a significant 
volume of model parameter space. 
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